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Abstract 

We study toy aging processes in hierarchically decomposed phase 
spaces where the equilibrium probability distributions are multifrac- 
tal. We found that the an auto-correlation function, survival-return 
probability, shows crossover behavior from a power law in the 
quasi-equilibrium regime {t <C t„) to another power law t''^ (A > x) 
in the off-equilibrium regime {t ^ tw) obeying a simple t/t^ scal- 
ing law. The exponents x and A are related with the so called mass 
exponents which characterize the multifractality. 
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1 Introduction 



The aging processes, i. e. relaxational processes to approach the thermal 
equihbrium, are extremely slow in glassy systems like spin-glasses and one 
can observe remarkable aging effects in experiments [|l], 0. Among the various 
phenomenological descriptions of the aging effects are due to Sibani and Hoff- 
man who proposed a scenario based on the concept, hierarchical diffusion 
[^]. The latter concept has been implemented in many toy models 
which we hereafter refer to as hierarchical diffusion models. The concept [|TU| 
is roughly as the following. Firstly, one considers that the free-energy land- 
scape consists of hierarchically nested valleys, which are usually described in 
terms of a certain tree structure. Then one introduces a relaxational dynam- 
ics in terms of a certain master equation which describes diffusion processes 
between different valleys driven by thermal hoppings over the barriers. Solv- 
ing the master equation, one obtains the time evolution of the distribution 
of probabilities to find the system at different bottoms of valleys (or leaves 
of the trees). 

In the present paper, we consider hierarchical diffusion in a class of trees 
which have the following two characteristics. Firstly, we consider that the 
backbone structures of the trees have self-similarity as in many of the previ- 
ously studied models. Secondly, we consider that the equilibrium probability 
distributions on the leaves are multifractal. The latter point is different from 
the previously studied models, whose equilibrium probability distributions 
are restricted to be uniform by their designs ( see however |^). We introduce 
the relaxational dynamics in terms of an exactly solvable master equation. 

We study aging processes after rapid temperature quenches in our model. 
An aging process appears as the growth of sub-trees in which the probabil- 
ity distributions are quasi-equilibrium (multifractal) whereas on larger scales 
than such sub-trees the probability distributions are still non-equilibrium 
(non- multifractal). As the result, a simple auto-correlation function, the 
survival-return probability, shows a characteristic crossover behavior: it de- 
cays by a power law t~^ in the quasi-equilibrium regime {t <C t^) but by 
another power law t~^ in the off-equilibrium regime {t ^ t„) and obeys a 
simple t/tw scaling. The exponents x and A turns out to be related with the 
so called mass exponents which characterize the multifractal properties of 
the probability distributions. 

The organization of this paper is the following. In section 2, we introduce 
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our hierarchical diffusion modeL In section 3, we analyze the aging effect in 
our model focusing on the scaling properties of an autocorrelation function. 
In section 4, we summarize this paper with some discussions. 

2 The Model 

2.1 Construction of Tree 

Let us consider a system of many states which have the following clustering 
property. Suppose that the system can be coarse-grained so that different 
states merge into fewer number of states. In Fig. 1, we show such a sys- 
tem represented as a tree. The magnitude of the resolution power increases 
downward along the h axis: the states (branches) are differentiated into more 
states (branches) as we raise the resolution power. We label the states dif- 
ferentiated with the maximum resolution power ( the leaves on the baselines 
of the trees) as Li (i — 1 .. . N) where N is the number of such states. 

For simplicity, we will only consider bifurcating trees: at every branch 
point we always have two branches stemming down. We will refer to the set 
of branches and branch points under a branch point, say A, as sub-tree A. 
We will refer one and the other side under a branch point A as side- / and // 
of A and present them on the left and right side respectively in the figures 
below. 

We construct the partition function of the system at the equilibrium as the 
following. The coarse-grained states (branches) at the same resolution level 
are considered to be all energetically degenerate with each other. However, 
the number of microscopic states they contain can vary, which means that 
the equilibrium local entropies and so the equilibrium statistical weights can 
vary as well. We consider that as a bifurcation takes place downward along 
the h axis (the direction to increase the resolution power), the equilibrium 
local entropy of a coarse-grained state is partitioned into those of two sub- 
states underneath. In order words, the partition function of sub-tree B is 
partitioned into those of side-/ and side-// sub-trees under B with a certain 
partition ratio, say 1 — ip{B) and ip{B) (0 < -0 < 1), respectively. 

Let us denote the child of B (a branch point just below B) on side-/ as 
Ci{B) and side-// as Cn{B). Then the ratio of the partition function of 
sub-tree Ci{B) and Cij{B) to that of B, which we denote as 7r(S, B), takes 
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the following values, 



(1) 



It is useful to generalize the above argument as the following. Suppose that 
a sub-tree D is enclosed by a larger sub-tree B. Let us denote the parent of 

D ( the branch point just above D) as Di and the grand-parent of D ( the 
parent of Di) as D2 and so on. Suppose that B is the i^''th ancestor of D, 
i. e. Dk = B. Then the ratio of the partition function of sub-tree D to that 
of sub-tree B, which we denote as 7r(L', 5), can be written as the product of 
tt's along the (unique) vertical path which connects D and B, 



For example, consider the set of leaves in a sub-tree B. Then 7r(Lj, B) asso- 
ciated with such a leaf Li can be interpreted as the relative statistical weight 
of the leaf among the set of leaves in sub-tree B. We denote the highest 
branch point as i?top and choose the partition function of the whole tree to 
be 1. Then the equilibrium statistical weight of a leaf Lj can be written as 
7r(Lj, i?top)- 

Here wc define some other terminologies for later uses. We denote the 
ancestors and descendants, which are the set of branch points above and 
below the branch point B, as A{B) and ^^{B) respectively. We denote the 
set of all branch points under side-/ and // of B as T>i(B) and T>ii{B) 
respectively (Note that Vi{B) U Vn{B) = V{B).) We denote the parent, 
grand-parent and the f^'th ancestor of B as Bi, B2 and B^ respectively as 
we already did above. We denote the lowest common ancestor of B and B, 
the branch point at the top of the smallest sub-tree which enclose both B 
and B, as A{B,B). For the convenience, we also introduce a hypothetical 
branch point -Bceiiing whose height is /iB^^in^g = 00 and set 7r(Stop, -BceiUng) = 1- 



We now introduce a stochastic dynamics of the temporal state in the hier- 
archically decomposed phase space. We consider that there is an thermally 
activated excitation associated with a branch point, say B, in the tree with 
which the temporal state can go from one to the other leaves in sub-tree B. 



7t{D, B) = 7t{D, Di)7l{D,, D2) . . . 7t{Dk-1, B). 



(2) 



2.2 The Master Equation 
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The excitations associated with the branch points at higher h are consid- 
ered to have higher activation energies to be excited. Thus we now redefine 
the vertical axis h as the scale of the activation energies of such excitations 
associated with the branch points. 

Wc now introduce a simple exactly solvable dynamics which describes the 
stochastic jumps between the leaves, i. c. the states which arc diffcrcntiable 
with the maximum resolution power. Here after we call the latter simply 
as states. Let us denote the probability to find the temporal state of the 
system at a state Li at time t as Pi{t). The time dependent distribution of 
the probability can be expressed in terms of a vector, 

p{t)^ip,{t),P2it),...,PN{t)), (3) 

which should become equal to the equilibrium probability distribution in the 
hmit i — > oo. We denote the latter as as p^^^ = {pT^tPI^, ■ ■ ■ where 

p^^ = 7r(L,,Stop). (4) 

Let us denote the transition probability to go from state Lj to Lj in a unit 
time as Wjj . Then the master equation for the evolution of the probabilities 
can be written as 

|p(i) = -Tp{t), (5) 

with 

-r,, = w,,-%i:Wfe,-. (6) 

k 

Note that the sum of the probabihty J2i^=iPi{t) is always conserved. 

We consider that the thermal jump process of the temporal state from a 

state (leaf) L to another state consists of two stages. In the first stage, the 
excitations associated with the branch points in A{L) (ancestors of L) are 
activated in a successive manner as the following. Suppose that the excitation 
associated with such a branch point B is activated. Then its parent Bi get 
a chance to become active or not with the probability exp(— (/i^^ — hs)) and 
1 — cxp(— (/?.B^ — hs)) respectively. If Bi becomes active, we repeat the same 
trial for B2. Otherwise, the successive ignition of the excitations stop there 
at the level hs- Thus the probability that the 1st stage ends at B is 

K 

w{B^L) = JJexp(-(/i^„-/i^„_J)(l-exp(-(/iB, -/is))) 

n=l 

= exp(-/iB) - exp(-/iBj (7) 



5 



where Aq = L and Ak = B. 

The 2nd stage is the falling down process from the height hs to a leaf 
of sub-tree B. Reminding that the leaves L G T>{B) have different relative 
statistical weights tt{L,B), we expect that the probabilities to fall into the 
leaves depend on their amount of local entropies in such a way that those with 
larger amount of local entropies have more chances to receive the temporal 
state. So we simply choose the probability to fall down to a state (leaf) L as 

w{L ^ B) =7i{L,B). (8) 

Combining the above two factors, we obtain the transition probability to go 
from L to L via B as 

w{L\B\L) = w{L^B)w{B^L) 

= 7r{L,B)[exp{-hB)-exp{-hB,)]. (9) 

Note however that such a process takes place only if both L and L belong to 
the sub-tree B. 

The transition probability Wij from Lj to Li is the sum of the transition 
probabilities over An of = 0, 1, . . . , M — 1 where Aq = A{Li, Lj) (the lowest 
common ancestor of Lj and Lj) and Am = -Bceiimg- Thus the off-diagonal- 
elements of the matrix W becomes 

A/-1 M-l 

^ij^j = w{Lj\An\Li) = J2 [exp(-^A„) - exp{-hA„+,)] 7r{Li,An), 

n=0 n=0 

(10) 

Using (^ and (PH]), it can be checked that the detailed balance condition 

w,,pf = w,,pr, (11) 

is satisfied with this choice. 



2.3 Solution of the Master Equation 

We now solve the master equation (^. The formal solution can be written 

as 

Pit) = exp(-rt)p(0) (12) 
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where p(0) is the initial distribution at t = 0. The probabihty that the 
temporal state which initially stay at Lj reach Li at time t is 



G.,(t) = [exp(-rt)],^.. (13) 

In order to calculate the propagator G, it is convenient to introduce a new 
matrix F with which we can rewrite F as 

F = (p<=^)^/2p(peq)-l/2^ (^4) 

where p'^'^ is the vector of equilibrium statistical weights p^'^ defined in (^). 
Note that F is a real symmetric matrix so that it has real eigen values. 

We now look for the N eigen states of F in a heuristic way. First of all, 
the static mode can be found easily as the following. The matrix F satisfies 
Fp<=i = due to (D and from which we obtain f (p^^)^/^ ^ g. The last 
equation means that the vector (p'^'^)^/^ is an eigenvector whose eigenvalue 
is i. e. static mode. There are N — 1 other eigen states {dynamic modes) 
left to be found. 

We construct here a set of vectors which consists of vectors localized under 
the branch points. We may call this set of vectors as umbrella set because 
of the localized shape of the amplitudes. On a branch point B we define a 
vector 

S,iB) = n'/'^{L„B)uiU,B) (15) 



where 



u{B,B) 




B 



B ^coiling 

B -Bcciling B eVi{B) 

B ^ -Bcciiing B e 'Du{B) 
otherwise. 



(16) 



Note that the vector S'i(-Bceiiing) is identical to the eigen vector of the static 
mode we obtained above. It can be easily checked that the vectors are nor- 
malized and orthogonal, 

Y.SmS^{B) = 6^^^. (17) 
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Since there are N — 1 eigen vectors on — 1 branch points and one static 
mode, which are all linearly independent with each other, they together con- 
stitute a complete set of dimension N. As shown in Appendix A, the vector 
Si{B) actually turn out to be the correct eigen vectors of T which have the 
associated eigen values 

z{B)=exp{-hB). (18) 
Note that ^(-Bcciiing) = (static mode) is endured since we have set /i_B^^iii„g = 

CXO. 

Now we rewrite some previously defined matrices in terms of the umbrella 
set. At first, T becomes 

t,, = j:S^{B)z{B)SAB). (19) 

B 

Then using the last equation and (|T5|) in (p!^), the matrix Tij becomes 

= Y.^iL„B)u{Li,B)u{L„B)z{B). (20) 

B 

Finally, we also rewrite the propagator G in terms of the umbrella set. Using 
(|13D, (|14D and (|l|) we obtain, 

G^jit) = {pTy^'Y.{S^{B)eM-mt)S,{B)}{pf)-'/' 

B 

= ^7r(L„5)u(Li,5)u(L,-,5)exp(-z(5)t). (21) 

B 

Using (^) and (|16[) and performing similar calculus shown in the appendix, 
we obtain the propagator in a more explicit form, 

M-l 

Gijit) = [exp(-z(A„+i)t) - exp{-z{An)t)] ^{U, An) 

n=0 

+ 6ijexp{-z{Ao)t) (22) 
where we defined Aq = A{Li, Lj) and Ai, A2, . . . , Am-i, Am = -Bcciimg- 
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2.4 Multifractality on Self-similar Trees 



2.4.1 Mass Exponents 

We consider in this paper self-similar trees on which the distributions of the 
equilibrium statistical weights have the following multifractal characteristics. 



Let us define the q-th moments of the statistical weights [0 as 



Mq{h)dh = J2 Kh- {hB-h^))'Ki{B,B)dh (23) 

BeV{B) 

where the over-line means the average over statistically independent sub-trees 
B. Suppose that a g-th moment have the following scaling behavior, 

Mg{h)dh ~ exp{T{q)h)dh, (24) 



where the exponent r(g) is called as a mass exponent [|T6l. If r(g) depends 
non-linearly on q, the distribution is regarded as multifractal. 

The geometrical self-similarity appears in the 0th moment Mo(/i), which 
is just the number of leaves of a tree of height h. The mass exponent r(0) 
is the fractal dimension of the tree and sometimes called as silhouette of 
the tree[^, which measures whether a given tree is slender or fat. A trivial 
exponent is r(l) which is always because of the normalization condition 
of the statistical weights. Another exponent which turns out to be quite 
important is the mass exponent of the 2nd moment r(2). As we see later, 
the two exponents r(0) and r(2) are related with the dynamic exponents of 
the survival-return probability in the present model. 



2.4.2 Randomly Branching Trees 

As an example of the trees which have the multifractal properties mentioned 
above, we construct here a specific class of random trees generated by the 
following randomly branching process (RBP). We use this specific example 
later when some demonstrations are necessary. 

An RBP starts with a single leaf, which is regarded as the highest branch 
point Btop and two branches stemming down from it. In a single step, the 
length of the branches under the lowest (new) branch points get longer by 
one unit, say db so that the height (from the baseline of the tree) of all branch 
points and so that the height of the tree get larger by db. The branching at 
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a leaf occurs with probability ^branch in a single step. When it takes place, 
the leaf becomes a new branch point and two new branches starts from it. 
Each of such events occurs independently from each other. Repeating this 
procedure, we obtain a backbone structure of a tree. 

Next we assign the weights on the tree. Consider a branch point B and its 
child (the branch points just below B) Cj on side-/ and Cu on side-//. We 
determine the variable ipB and assign n{Cj, B) = 1 —ips and 7r(C//, B) = ipB 
to the branches on the side-/ and // of B. The variable ipB on each branch 
point is chosen randomly from the distribution 

F{ip)dil) = probability that ipB lies between ip and ip + (25) 

We perform this procedure for the whole branch points. 

Consider a tree which grows larger by the RBP. A natural consequence of 
the RBP is that the backbone structure possesses statistical self-similarities. 
On the other hand, the statistical weights on the leaves are successively par- 
titioned further into more and more fine pieces by the RBP. It is well known 
that if such a process is repeated, one oftenly finds very peculiar distribution 
of the weights : some set of pieces which have rather larger weights but neg- 
ligibly smaller population compared with the typical ones come to rapidly 
dominate the total sum of the statistical weights (the partition function of 
the whole tree) as the branching proceed further. This phenomena is called 
as curdling [|T5| of multifractal objects. 

As we show in Appendix B, we actually obtain the scaling property of 
the form (|2^) in the case of the random trees generated by the RBP. The 
mass exponent is. 



It can be seen that it is generally non-linear with q. In Fig. 2 we show an ex- 
ample of r(g) on random trees generated by a RBP. A special case when T{q) 
becomes linear with q is when the following two conditions hold: Pbranch = 1 
(deterministic branching) and F{iIj) = 6{ijj — 1/2) (always symmetric parti- 
tion) . 

The real samples of such random trees can be generated numerically by 
by the following Monte Carlo Method. In one Monte Carlo step (MCS), 
the height of all branch points are raised by db. Simultaneously a pseudo 
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dijFiij) {r + (1 - + (1 - Pbranch) 



(26) 
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random number is generated for every leaf and if it is smaller than pbranch, a 
bifurcation takes place: the leaf becomes a new branch point and two new 
branches start from it. Fig. 1 is actually an example of such a random tree 
obtained by simulating the RBP of Pbranch = 0.10, db = 1.0 and F^ip) = 
5{ip - 0.2) for 20 MCS. 

3 Aging Effect 

3.1 Growth of Quasi-Equilibrium Domain 

We now consider an aging process after rapid temperature quench from high 
temperature. For this purpose, we choose the initial configuration as, 

P(0) = ^. (27) 

After waiting for (waiting time), the probability distribution evolves 
up to, 

p(tw) = G(t„)p(0), (28) 

which eventually become p*^"^ a.st„ ^ oo. Note that the initial non-equilibrium 
distribution is not multifractal because it is uniform, while the final fully- 
equilibrated distribution is multifractal. Hence the aging process in the 
present context can be understood as the process to approach a multifractal 
distribution from a non-mult if r act al distribution. 

In order to see how the system age, it is convenient to define 

ri(t„) =p,(tw)/pr= ^Eexp(-5(5)t^)7r-i(5,5top)M(L„5), (29) 

where we used (^) in ( pSD and defined 

u{U,B)=u{U,B) u{Lj,B). (30) 

jeV{B) 

In Fig. 3 we show an example of ri{t^) calculated using the exact solution 
of the master equation solved on a real sample of random tree shown in 
Fig. 1. One can see that, as increases, Viit^) of different states come to 
join with each other successively and constitute groups, among each of which 
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the values of ri(t„) are common. Note that as far as the transitions within 
such groups are concerned, the detailed balance (|lT]) is fulfilled. Thus we 
may call such a group of states as a quasi- equilibrium domain. 

In order to understand the growth mechanism of the quasi- equilibrium 
domain in a more formal way, let us introduce a characteristic height h^^{t„) 
which grows logarithmically with t^, 

h'^'it^) = log(t„). (31) 

Due to the factor exp{—z{B)t^) ~ exp(— exp(/i°^(tw) — hs)) , the contribu- 
tions from the branch points lower than h'^^{t^) in the r.h.s. of the equation 
( p9D are negligibly small, compared with those from the branch points higher 
than /i"^(tw)- So it can be roughly approximated as 

r,(t„)^ expi-z{B)t^)n-\B,Btop)u{U,B). (32) 

/i(B)>/i<=ff(iw) 

Consider a pair of states Lj and Lj whose lowest common ancestor is A{Li, Lj). 
Suppose that after certain waiting time t^, the characteristic height h'^^itl^) 
becomes larger than hA{Li,L )■ Then from the definitions (|30|) and ([T6|) , all 
the terms that survive in the sum of ( p2D become the same for Lj and Lj. 
Consequently Viit^) = rjit^) holds forever for 3> t^. 

To summarize, the aging process of the present hierarchical model is un- 
derstood as the growth of aged sub-trees or quasi-equilibrium domains. Here 
we mean by a quasi-equilibrium domain a sub-tree under a branch point lower 
than the characteristic height h'^^lt^), which grows logarithmically with t^. 
The probability distributions inside an aged sub-tree is almost the same 
as that of fully-equilibrated one except for a common multiplicative factor. 
This is one of the most important consequence of hierarchical diffusion and 
actually found to be the case in the relaxational dynamics of microscopic 
spin-glass models |]lT] . 



The probability distribution inside an aged sub-tree is multifractal and 
possess the same scaling behaviors as that of the fully-equilibrated one while 
on higher scale h ^ h^^{t^), the system is still highly non-equilibrium or 
young in the sense that the probability distribution is not multifractal. 
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3.2 Survival-Return Probability 

We now introduce one of the simplest auto-correlation function as a prove, 
that is the survival-return probability, 

q{t + t,, t„) = ^ pi{t^)G{t)ii. (33) 

i 

It measures the probability that the system returns to the leaf where it stayed 
at time tw after additional traveling of t. 

Using (pTD, the autocorrelation function can be rewritten in terms of the 
umbrella set as 

q{t^,t + t^) = ^pi(tw)G(t)ii 

i 

= ^exp(-ZBt) ^pi{t^)-n{Li,B)u^{Li,B) 

B i 

= J dzpt^{z)exp{-zt), (34) 
where we defined a kernel pt^ (z) as 

ptjz) = J2^i^- ^b) |5:Pi(tw)7r(L„i?)n2(L„i?)| . (35) 

In the latter subsections, we study the scaling behaviors of the survival-return 
probability in random trees generated by RBP focusing on the role played 
by the waiting time t^. 

3.3 Two Extreme Cases 
3.3.1 Zero Waiting Time 

We consider at first a special case of zero- waiting time case t„ = 0, which 
means that the system is in an extremely non-equilibrium condition at t = 0. 
Since we have set the initial condition as (0), we obtain 

B 

= ^T.Hh-hB))dh, (36) 

B 



13 



where we used (0) and defined a variable h = — log(2;). 
Using dH) and (H), we obtain 

LY^6{h-hs)dh = ^ E 6{{hs,^^-h)-{hs,^^-hs))n%B,B,,,)dh 

B BeV{Btop) 

~ j^Mo{hB,^^-h)dhc:^exp{-T{0)h)dh 

~ z'-^dz, (37) 

where s = r(0) is the silhouette. In the last equation we used ~ MQ^hs^^p) 
where -Stop is the highest branch point. In the above equations, we approx- 
imated the sums by their mean values (^) assuming that the contributions 
of the deviations from this mean value vanish in the thermodynamics limit 
N — »• cxo, i.e. self- averaging. This assumption is valid on trees generated by 
the RBP because quantities on sub-trees under different branch points at the 
same height are statistically independent from each other. 

Then using (0), we obtain a power law decay in the off-equilibrium limit 

g(0,t) = J dzn{z)exp{-zt) (38) 

where the exponent A is equal to the fractal dimension or silhouette of the 
tree, 

X = s = r(0). (39) 

The above result is similar to those of the previously studied hierarchical 
diffusion models, which also yield power law decays whose exponents are 
related with the silhouette of the trees It is, however, not surprising 

because the distributions of the equilibrium probabilities in such models are 
uniform which is also the case for the present choice of the initial condition 



3.3.2 Infinite Waiting Time 

The special case of infinite waiting time t„ = oo is also of interest. In this 
limit, the system is fully-equilibrated or aged at t = 0. Using Pi(oo) = p'^^ 
and (H), we obtain 

Pooiz) = Y.Si^-m)\T.PiML^,B)u\L,,B).] (40) 
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^ 5{z - z{B))7i{B, 5top) E B)u\L,, B) 



B 



Y,5{z - z{B)y{B,B,,^)^B{z), (41) 



B 

where we defined 



= 5{z-zb) E 7,\L,,B)u\L,,B) (42) 

L,(LV{B) 

'ipB 



1 - tp^ 



-6{z-zb) E ^'(^.,5) 

^ L,eVi{B) 



+ E (43) 

It is sufficient to consider the scaling property of the 1st term in the last 
equation. We can rewrite it as, 

5iz - zb) E ^'(^- B) ~ M^{h) ~ z-<^\ (44) 

L,eVi{B) 

where we wrote h = — log(2;). In the ffist equation, we evaluated the sum by 
the mean value ( ^^1) assuming self-averaging property and used (^^. Thus 
^b{z) scales with z as $^(-2) ~ z~'^^'^\ In the same way we obtain, 

E'5(^ - zb)'k{B, Btop)dz ~ MiQiB,,^ - h)dh 

B 

dz , . 

~ const — (45) 

Z 

where we wrote again h = — log(2;) and used r(l) = 0. Combining above 
results we obtain 

Pooiz)dz ~ z-^^^'^-^dz. (46) 

Using ( |4^ ) in (Q), we finally obtain another power law in the fully-equilibrated 
limit, 

q{t + 00, 00) ~ t"^, (47) 
where x is an exponent defined as, 

X = -r(2). (48) 
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Let us make some comments on the difference of tlie two dynamical expo- 
nents A and x. Generally, the inequality \ > x and equivalently r(0) > — t(2) 
hold as long as r(g) decreases monotonically with increasing q and concave 
downward. (Note that r(l) = holds always as we mentioned before.) The 
latter conditions seem to be usually satisfied. This is certainly the case on 
the mass exponents of trees generated by the RBP whose formula is given in 
(|26D . The equality holds only when Pbranch = 1-0 (deterministic branching) 
and F{il)) = 6{ip — 0.5) (always symmetric partition) so that the equilib- 
rium probability distribution becomes the same as the initial non-equilibrium 
(non-multifractal) distribution. 



3.4 Crossover Behavior : Aging Effect 

Now we consider the case of finite waiting time, in which we expect some 
waiting time effects, i. e crossover from quasi-equilibrium to off-equilibrium 
behavior. It is now convenient to introduce another kernel p{z, z) such that 



Ptw(^) = / dzp{z,z)exp{-zt^). 



(49) 



Then the autocorrelation function f|33) can be rewritten as 



q{t + t^,t^) = J dz J dzp{z,z) exp(— 2;t) exp(— it„). 



(50) 



From (^), (^) and (^T]) we obtain the explicit form of the kernel p(z, z) as 



B B 



7r(L„ B)n{L,, B)u\L„ B)u{U, B) ^ u{U, B) \ . (51) 



The scaling form of p(z, z) is studied in Appendix C. Here we read the result, 

, dz dz 



p{z, z)dzdz 



z z 

zY^^^i , ,^^dzdz , 

z r < z). 

zj z z 



(52) 
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Using ( p2D in (|50|), we finally obtain 

q{t + U, U) ~ r^'^i{t/U) + t-^(°)t;(°)+^(')g2(t/tw), (53) 
where we defined 

/ny/ s 
dyy'^^^^'^ exp{-y) J dyy'^ exp{-y) 

Us) = Cj dyy^^'^~' exp{-y) f dyy-^^^'^+^^'^'>-' exp{-y) , (54) 

J Jy/s 

where Ci and C2 are numerical prefactors. 

From the above results, we find that the autocorrelation function obeys 
the following simple t/t^ type scaling, 

g(t + t„,tw) ~t~"g(tAw), (55) 
where the scaling function q{y) behaves as 

with X = — r(2) and A = t(0). 

From the above scaling from, it can be seen that the autocorrelation func- 
tion crossovers from quasi-equilibrium behavior to off-equilibrium behav- 
ior t^^ at around t ~ t„. This crossover behavior appears due to the growth 
of the quasi-equilibrium domain in which the probability distribution is mul- 
tifractal while on larger scale, it is still non-equilibrium (non-multifractal). 
The inequality of the two exponents X > x means that the off-equilibrium 
decay is faster than the quasi-equilibrium decay, which is intuitively satis- 
factory. 

We show in Fig. 4 some examples of the crossover behavior of the autocor- 
relation function, which was obtained using the exact solutions of the master 
equation on real samples of random trees. The random trees are generated by 
Monte Carlo method which simulate the RBP of a) Pbranch = 0.10, db = 1.0 
and Fitfj) = 5(^-0.2) and b) planch = 0.10, db = 1.0 and F{ij) = 6{ij-0.5). 
The random average was took over 10^ different realizations of such trees gen- 
erated by 50 Monte Carlo steps. The predicted power law and with 
a) A = r(0) = 0.0953... and x = -r(2) = 0.0325... (see Fig. 2) and b) 
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A = r(0) = 0.0953... and x = — t(2) = 0.0513... which are obtained from 
(26), are also included in the figure. The curvatures of the curves at lower 
values of q are due to the finite size effects. In Fig. 5 we show the scaling 
plot of the data shown in Fig. 4 a). The curves of different are plotted 
against and shifted vertically so as to converge to a master curve. One 
can well see that the data are consistent with the predicted scaling laws (^) 
and (^). 



4 Discussions 

We have studied aging effects in a simple exactly solvable model of hierarchi- 
cal diffusion. We considered the case that equilibrium probability distribu- 
tion has multifractality. A specific way to generate such trees by randomly 
branching processes (RBP) are introduced for demonstrations. Aging pro- 
cesses after temperature quenches appear as the growth of aged sub-trees in 
which the probability distribution is in quasi-equilibrium and multifractal. 
In the thermodynamics limit, the height of the tree becomes infinite and the 
true equilibrium cannot be attained in any large but finite waiting time t^. 
Consequently, the waiting time dependence persists for the whole range of 

except for t„ = cxd, i. e the ergodicity is weakly broken |jl2[. We found 
that these properties clearly reflect in the survival-return probability and 
brings about the characteristic crossover from quasi-equilibrium behavior to 
off-equilibrium behavior. 

Let us make some comments on the robustness of the scaling properties 
of the autocorrelation function we obtained in our exactly solvable model. 
Note that there are other possible choices of the transition matrix other 
than our present choice, which describe hierarchical diffusions and endures 
the detailed balance condition (pT]). For instance, one may define another 
transition matrix by replacing nlLi, An) in (jl^) by 7r(Lj, An)~^. One can also 
construct transition matrices considering that transitions between a pair of 
leaves occur only over their lowest common ancestor. We investigated the 
solutions of the master equations with these alternative transition matrices on 
random trees generated by the RBP by numerical diagonalizing the transition 
matricies. Interestingly enough, we found that the scaling behaviors of the 
autocorrelation function appears essentially the same as that of the exactly 
solvable one presented in this paper and one only needs some renormalization 
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of the global unit of time. This fact implies that the scaling properties are 
robust to some extent against minor changes of the model. 

It is interesting to note that the crossover behavior from quasi- equilibrium 
behavior (t <^ t„) to off- equilibrium behavior {t ^ t„) obtained in the preset 
toy model, is very similar to that observed in the relaxational dynamics of 
some micro-scopic models of random systems such as 3D spin-glass model |]13| 



and 1+1 dimensional directed polymer in random media |Tj]. In the latter 
models, the autocorrelation functions obey t/t^ scaling law with two power 
law decays at the quasi-equilibrium regime and at the off-equilibrium 
regime, which is what we found in the present phenomenological toy model. 
Thus it is tempting to speculate that our present toy model will provide 
a clue to understand the link between phenomenological pictures based on 
hierarchical diffusion and the glassy dynamics of realistic systems [|18|. 
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A The Umbrella Set 



In this appendix we show that the umbrella set defined in (p!5[) , ([T6|) and (|T^) 
are the true eigen states of the matrix F defined in (j^). We prove this by 
checking if the umbrella set correctly reproduce the original transition matrix 
defined in (|10D. 

Using (|T5|) and (|18]) in (|20|) we obtain 

-Ti^j = -Y.n{U,B)uiLi,B)z{B)uiLj,B) 

B 

= 7T{Li, A{Li, Lj)) exp(-/iA(L„L,)) 

J2 7i{L„B)exp{-hB)u\L„B) 

BeA{A{L„Lj)) 

M 

= 7r(Li, Ao) exp{-hAo) - ^(^i' exp(-/iA„)M^(^i, ^n) 

n=l 

where we defined Aq = A{Li,Lj) and Ai, A2, . . . , Am-2, ^m-i = B^op and 

Am = -Bceiling- 

We further rewrite the right hand side of the last equation as the follow- 
ing. The factor u{Lj, An) in the last equation can be replaced by u{An-i, An) 
due to the definition (p!6|). And the factor 7r(Lj,A„) can be decomposed as 
7r(Lj, An~i)7i{An-i, An). Then we can use the identity 7r(A„_i, An)u'^{An_i, An) 
1 — TT{An-i, An), which follows from (16) and (0). Then we obtain 

r. h. s. = 7i{Li,Ao)exp{-hAo) 

M 

- 7i{Li, A„_i) (1 - 7r(y4„_i, An)) exp{-hA„) 

n=l 
M-1 

= ^ [exp(-/iAj - exp(-/iA„+j] vr(Li, A„) (57) 

n=0 

where we used = ^Sccmn = cxd in the last equation. Then using the 
relation (^, we see that the off-diagonal elements of the transition probability 
Wj^j defined in equation ([10|) is correctly reproduced by the umbrella set. 

Next we check if the matrix F written in terms of the umbrella set, prop- 
erly conserve the total probability. Taking the sum over i of both sides of 
we obtain 

Er., = Y(T.^iL^,B)u{L,,B)]z{B)n{L„B) 

i B \ i J 
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= (58) 

The last equation follows from the following. For B = -BcciUng, the contri- 
bution is since i(-Bcciiing) = 0. And for B ^ -Bceiiing, one finds again zero 
contributions using (p!5D, 



(1 - ^Pb) - \ ■ = 0. 



V 1 - V^B \ i^B 

Thus we obtain the last equation of (P3). 



B Mass Exponents on Random Trees 

In this appendix we study the scaling property of the g-th moment of the 
probability distribution on the random trees generated by the randomly 



branching process (RBP) and derive the formula (pBI) . 

Let us denote the probability that the q-th moment Mq{h) of a random 
tree of size h = mdh takes value x as ujq{m, x). Considering that larger trees 
can be constructed by smaller sub-trees, we obtain the following recursion 
relation for u!q{m,x), 

ujq{m + 1 

5 'f^) Pbranch J dyidy2 d'ijjF{ip)LUq{m,yi)ujq{m,y2) 

X (5(V'V + (1 - V^)^l/2 - X) + (1 - Phr^nch)uJq{m,x). (59) 

In order to solve this integral equation, it is convenient to introduce a 
generating function defined as 

Zq(u,m) = j dxex^{ux)ujq{m,x). (60) 

The expectation value of the q-the moment can be obtained as 

d 



Mq{h) ~< X >q^m= J dxXUq{m,x) = —Zq{u,m] 



(61) 

u=Q 
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Multiplying exp(-ux) on both sides of ( pQ]) and integrating over x, we obtain 
the recursion relation for Zq{u, M) 



Zq{u,m + 1) = ^branch / dljj F (ip) Zg{u'ljj'^ , TTl) Zq{u{l - l/j^ , TTl) 

Jo 

+ {1 - PhrAnch)Zq{u,m). 

Then we obtain the recursion relation for < a; > 



(62) 



q,mi 



< X >q,M+l 



Pbranch Pbranch / 

"'0 



Solving the last equation with < x >g,i= 1, we get 

Mq{h) ~< X >q,M= exp(r(g)/i), 
where the mass exponent r(g) is obtained as 



< X >q,M ■ 

(63) 
(64) 



r(g) = db ^ log 



1 

Pbranch 

/ #F(V') {V-^ + (1 - ^)«} + (1 Pbranch J 





(65) 



In order to further investigate the multifractal properties, it is convenient 
to introduce the exponent of singularity a defined as 



7c{B,B)=exp -a{B,B){hB - h 



B) 



(66) 



Then one can obtain the distribution of a or /(a) spectrum using the well 



known procedure [16| and discuss curdling. However we don't discuss about 
it here [|l^ . 



C Scaling form of Kernel 

In this appendix we study the scaling property of the kernel p{z, z) with 



respect to z and z. Its explicit form (pll) is 



p{z,z) 



lEE5(.-5(5))5(i-5(i?)) 



N 



B B 



J2 AU. B)n{L,, B)u\U, B)u{L„ B) ^ n(%, B) . (67) 
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Note that, from the definitions (|I6|), the factor v?{Li, B)u{Li, B) is non-zero 
only when the leaf Li is under both B and B. 

At first we consider the case z > z. In this case, the terms which survive 
in the sum (|67|) are those in which B is an ancestor of B and B is an ancestor 
of Lj. Let us introduce h = — log(2;) and h = — log(i). Then we obtain 

p{z,z)dzdz {z>z) = J2 S{z-z{B))n\B,B,,p) 

B&V{Btop) 

X Yl Kz-z{B))^\L„B)u{L,,B) 

Lj&ViB) 

X S{z-z{B))6{z-z{B))7r\B,B)u{B,B) 

B£V{B) 

X Y S{z- ziB))TT^{Li,B)u^{L„B)dzdz 

~ ^Mo(/iB,„^ - h)Mo{h)M^{h - h)M2{h)dhdh 

~ ,-r(2)^^ (68) 
z z 

where we evaluated the sums by their mean values (^), assuming self- 
averaging property, and used p6|) and N ~ Mo^hB^^p)- 

The other case z < z can be analyzed in the same way. Considering that 
the terms which survive in the sum ( |67D are those in which B is an ancestor 
of B and B is an ancestor of Li, we obtain 

p{z,z)dzdz {z<z) = ^ Y S{z - z{B))7T°{B,Btop) 

Y 5{z - z{B))S{z - z{B))7i\B, B)u''{B, B) 

BeV{B) 

Y 6{z-z{B))n\L,,B)u{L,,B) 

Y 6{z - z{B))'K^{Li,B)u{L,,B)dzdz 

Lj€V{B) 
1 

N 

z 



X 



X 



X 



:Mo(/iij,„p - h)Mi{h - h)M2ih)Mo{h)dhdh 



(69) 

z z 
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where we used again N ~ Mo(/i_Bt^p). Combining above results, we obtain 
31. 
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FIGURE CAPTIONS 



Fig. 1 Hierarchical organization of states. This example is generated by a 
RBP explained in section |2A^ (pbranch = 0.10, db = 1.0, F(V') = - 0.2), 



20 MCS.) All the branches on the left side have weight 1 — ip = 0.8 while 
those on the right side have weight ip = 0.2 as indicated in the figure. 

Fig. 2 The non-linear behavior of the mass exponent r(g) v.s. q: the curve is 
obtained using (|6|) for RBP of pbranch = 0.10, db = 1.0 and F{tp) = 5(^-0.2). 
The two important values of r(g) at g = and q = 2 are indicated by the 
arrows for later reference. 

Fig. 3 Growth of local equilibrium domain with increasing waiting time t„ : 
the plot of rj(t„) = Pi(t^)/Pi^ vs. t^ on the random tree shown in Fig. 2. 

Fig. 4 The crossover behavior of q{t^ + 1, t^) at different waiting times t^ on 
two specific kinds of random trees a) and b) (see text) . The waiting time t^ 
varies as t^ = 0, 10, 10^, . . . , 10^^, oo from the lowest to the top curve. 

Fig. 5 The t/t^ scaling plot of q(t^+t, t„) : the data of of different t^ presented 
in Fig. 5 a) except t^ = and oo, are used in this plot. The vertical scale is 
arbitrary. 
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¥=0.2 





Fig. 5 



LIST OF CHANGES MADE IN THIS REVISED VERSION 



The subsection which was titled 'Distribution of Relaxation Times' in 

was 



the original version, which was placed just after subsection p^ , 
removed and the necessary part of it is now included in subsection 



The figure 1 of the original version, which was intended to explain the 
idea of hierarchical decomposition of the phase space, was removed. 
This is because figure 2 in the original version (figure 1 of the revised 
version), which shows an example of random trees generated by RBP, 
can be used to explain the idea. 

The list of references were updated and some additional references ( 
the last three of it) are added. 

Minor changes such as corrections of miss typing and refinements of 
some sentences are made. 



